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Abstract 

00 

£N| Let P = P(h) be a semiclassical pseudodifferential operator on a Rie- 

mannian manifold M. Suppose that u(h) is a localised, L 2 normalised 
' family of functions such that P(h)u(h) is O(h) in L 2 , as h — > 0. Then, 

'"""j for any submanifold Y C M, we obtain estimates on the LP norm of u(h) 

"^-n restricted to Y, with exponents that are sharp for h — > 0. These results 

pH generalise those of Burq, Gerard and Tzvetkov [4] on LP norms for restric- 

tion of Laplacian eigenfunctions. As part of the technical development we 
prove some extensions of the abstract Strichartz estimates of Keel and 
S Tao 0. 

Let P = P(h) be a semiclassical pseudodifferential operator on a Rieman- 
nian manifold M. We will assume that P has a real principal symbol, and that 
( ~^ ) its full symbol is smooth in the semiclassical parameter h. Other more technical 

assumptions on P are given in Definition |1.6| We prove estimates for approxi- 
mate solutions u — u(h) to the equation P(h)u(h) = 0. As usual in semiclassical 
O^l analysis we assume that u{h) is defined at least for a sequence h n tending to 

zero. 

Our precise definition of approximate solution, or quasimode, is that P(u) = 
< 3\ Oli(K) as h — ► 0. This definition is natural with respect to localisation: if 

P(u) = Oui(K), and x is a pseudodifferential operator of order zero (with a 
symbol smooth in h), then P(x u ) is a l so 0^2 (h). We will make the assumption 
that u{h) can be localised, see Definition [T3] and therefore will be able to reduce 
the problem to one of local analysis. 

Given a submanifold Y of M, we estimate the LP norm of the restriction of 
u to Y, assuming the normalisation condition ||u||l 2 (a/) = 1- These estimates 
are of the form < Chr s where 8 depends on the dimension n of M, 

the dimension k of Y and p (except for one case where there is a logarithmic 
divergence) — see Theorem |1.7| In every case the exponent S(n,k,p) given by 
Theorem 1 1.7| is optimal. Figure [T] shows the exponent S for a hypersurface and, 
for comparison, the LP estimates over the whole manifold (Sogge llj for spectral 
clusters and Koch-Tataru-Zworski [S] for semiclassical operators). Figure [2] 
shows S(n,k,p) for submanifolds of codimension greater than one. 
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Figure 1: S(n, k,p) plotted against 1/p for the hypersurface case (solid line) and 
whole manifold estimates (dashed line) 
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Figure 2: S(n, k,p) plotted against 1/p for submanifolds of codimension greater 
than one 
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The potential growth/concentration of the quasimodes of a semiclassical 
operator is of great interest due to the connection to Quantum Mechanics. It is 
from Quantum Mechanics that we get the important set of motivating examples, 

P = h 2 A g + V (x) (1) 

here A g is the (positive) Laplacc-Bcltrami operator associated with the metric 
g. We can transition between this picture and the usual eigenfunction picture 
of Quantum Mechanics by dividing the eigenfunction equation 

Am + V (x)u = Eu 

by E. Then setting E = l/h 2 we have 

h 2 Au + h 2 V (x)u -u = 

or 

Pu = 

where P is as in with a potential term of h 2 Vo(x) — 1. Therefore the higher 
eigenvalue asymptotics of eigcnfunctions of Quantum Mechanical systems cor- 
responds to the h — > limit in semiclassical analysis. When Vq(x) = this 
problem reduces to estimating the size of Laplacian eigenfunctions restricted 
to a submanifold. A complete set of estimates for Laplacian eigenfunctions on 
compact manifolds is given by Burq, Gerard and Tzvetkov [4]. 

A number of different techniques for studying the potential concentrations of 
eigenfunctions are available. A large body of recent work focuses on semiclassi- 
cal measures (see for example Anatharaman [I], Gerard-Leichtnam [B], Zelditch 
[T2"] and Zelditch-Zworski [13]) ■ Sogge's work [TT] on spectral clusters give es- 
timates for |H| iP ( M ) of the form |M| LP ( M ) < \~ s ( n 'P) where A is the eigenvalue 
of u. This work is extended into the semiclassical regime by Burq, Gerard and 
Tzvetkov (3j (for the Laplacian) and Koch, Tataru and Zworski [9] (for semi- 
classical operators) . Multilinear estimates for spherical harmoics have also been 
obtained by Burq, Gerard and Tzvetkov [5] In related work Koch and Tataru [5] 
give LP estimates for eigenfunctions of the Hermite operator H = —A + x 2 . In 
2004 Reznikov [TO] proved bounds for restrictions of Laplacian eigenfunctions to 
curves where the underlying manifold was a hyperbolic surface. In 2007 Burq, 
Gerard and Tzvetkov [1] produced results giving LP estimates of the restriction 
of eigenfunctions of the Laplace-Beltrami operator on a compact manifold to a 
submanifold. This work directly extends these results using techniques found in 
Koch-Tataru-Zworski [9 and Burq-Gerard-Tzvetkov [2] to move them into the 
more general semiclassical setting. 

To continue we must define some objects from Semiclassical Analysis and 
give some basic results. A more detailed discussion of Semiclassical Analysis 
can be found in [3], [5] and [2], however for the reader's convenience the main 
definitions and results used in this paper are provided in Section [l] 
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1 Semiclassical Analysis 

Scmiclassical analysis allows us to study Pseudodiffcrcntial and Fourier Integral 
Operators depending on a parameter which we denote as h. We think of this 
parameter as being small and obtain error terms bounded by powers of h. As 
in the normal pseudodiffcrcntial calculus an operator P acting on L? functions 
u is given by its symbol p(x, £, h) and a quantisation procedure. 

Definition 1.1. Let p(x,£,h) <E C°°(T*R n ) be a symbol in the symbol space 
S m . We define the left semiclassical quantisation p(x , hD) as 



Remark 1.2. For real symbols the Weyl quantisation p w (x , hD) is self- adjoint. 
For this reason it will sometimes be more convenient to use p w (x, hD) in place 



Definition 1.3. A function u depending parametrically on h is said to satisfy 
the localisation condition if there exists \ € C C (T*M) such that 



where S is the space of Schwartz functions. 

This assumption allows us to move from a global problem to a local one. As 
X has compact support in T*M we can write 




and the Weyl semiclassical quantisation p w (x , hD) as 




ofp(x,hD). 



u = x(x,hD)u + O s (h oc ). 



N 
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for some N < oo where each \i has arbitarily small support. As noted previ- 
ously the notion of an approximate solution is preserved under such a localisa- 
tion. Now we may assume that we are working on a coordinate patch of T*M. 
Therefore we identify M with R n and Y with R k . An element x G M will be 
denoted x = (y, z) where Y = {z = 0}. An element £ £ T*M will be written as 
£ = (£y,£z)- Note that if M is a compact manifold the localisation requirement 
in the spatial variables is trivially satisfied. 

As we assume that p(x, £, h) is smooth in h we can write p{x, £, ft-) = Po(£) £)+ 
hq(x,£,h). Now as u is localised, 

P(x,hD) X {x,hD)u = L 2(h) => P (x 1 hD)x(x,hD)u^ L *{h). 

For the rest of this paper we will therefore assume that we are working with a 
symbol p(x,£) independent of h. 

Using the localisation assumption we are able to get a bound on \\u\\ LP in 
terms of \\u\\ Lq where q < p. We have 

u = x(x,/iL>)u + 0(/i 00 ) 

= h- n j K (x, ^) u(y)dy + 0(h°°) 

where 

K(x, Z ) = J ^J X (x,0e i <^>^. 

A bound of c(l — \x— z\)~ N is found for \K(x, z) \ by repeated integration by parts 
and by then applying Young's inequality the following estimate is obtained. 

Lemma 1.4. Suppose that a family u = u(h) satisfies the localisation condition 
then for 1 < q < p < oo 

Ml* < h n{1/p - 1/q) \u\ Lt + 0{h°°). 

In a couple of places we will want to use this estimate over a submanifold 
rather than the full manifold. To do this we require localisation to hold if some 
variables are fixed. 

Lemma 1.5. If u satisfies the localisation conditions then there exists some 
x(y, £, v ) compactly supported such that 

R Y u = X (y, hD y )R Y u + 0(h°°) 

where Ry is the restriction operator onto the submanifold Y. 

Proof. First as u(x) is localised we can replace u(x) with x( x , hD)u. Let 
ip(y, £y) e C^(M fe xM fe ) such that ip(y, £„) = 1 for all (y, £ y ) such that (y, z, £ y , £, z ) € 
Suppx(a;,£) for some (z,£ z ). As = ip(y,£, y )x( x >0 repeated application 

of non-stationary phase gives 
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(Id - ^{y, hD y ))R YX (x, hD)u = 0(h°°) 

which gives 

Ryx(x, hD)u = i/)(y, hD y )R Y x{x, hD)u + 0(h°°) 

as required. 



□ 



Using this localisation condition we can prove that when p{x, £) is bounded 
away from zero the local contribution is small. From (|9J, Lemma 2.1) we have 
that if > 1/C on a local patch then we can invert p(x,hD) up to 

order h°° . That is, choosing x( x >0 supported on this patch, we can find some 
q(x, hD) such that 

q{x, hD)p(x, hD)x{x, hD) = X (x, hD) + L 2^ L 2 (h°°) 

and 

p(x, hD)q(x, hD) X (x, hD) = X (x, hD) + L 2^ L 2(h°°). 
So if p[x, hD)u = Otfih) and |p(x, £)| > 1/C we can invert p(x,hD) to get 

X (x,hD)u=0 L2 (h). 



Now using Lemma 1.4 to estimate ||u|| Loe by \u\ L a we have 

ML = o(^ ( - 2)/2 ). (2) 



To get the L 2 norm of the restriction of u to Y we use Lemma 1.4 again this 
time only in the z coordinates. We have 

My, 0)|| i? < \\u(y, z)\\ LrL2y <h-^ \\u(y, z)| l;l , . (3) 

So the L 2 norm of u when restricted to a submanifold is 0(h 2 ). Inter- 
polating between Q and ([3| gives us better L p estimates than those given by 
Theorem 1.7 Consequently we can ignore regions where p(x, £) is bounded away 
from zero. 

This reduces our problem to localising around points (xo, £0) where p(xq, £q) = 
0. To proceed we need to place some non-degeneracy conditions on p(x,£). 

Definition 1.6. A symbol p(x , £) is admissible if it satisfies the following non- 
degeneracy conditions: 

(Al) for any pair (xo,£o) such that p(xo,£o) = 0; d£P(xo,£o) 7^ 

(A2) the second fundamental form on {£ | p(xo,£) — 0} C T* Q M is positive 
definite. 
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The first condition will be used to convert this problem into one regarding 
evolution operators. The second condition is needed for some later stationary 
phase estimates. The main result of this paper is below. 

Theorem 1.7. Let (M,g) be a smooth Riemannian manifold with no boundary 
and let Y be a smooth embedded submanifold with dimension k. Let u(h) be a 
family of L? normalised functions that satisfy Pu = 0^2 (h) for P a semiclassical 
operator with symbol p(a;,£). Assume further that u satisfies the localisation 
property and that the symbol p{x , £) is admissible. Then the LP norms restricted 
to Y are: 




2n < n < oo 

5(n,n-l,p) = { \ V n - 1 - P - (4) 

and for k ^ n — 1 

*(n,fc,j>) = {==i-*,2<p<oo, (k,p)jL(n- 2,2). (5) 

For k — n — 2 the L? estimate is 

\u\ L , {Y) <h~^\oil 2 (\lh). 

Remark 1.8. Apart from the log loss in the (k,p) = [n — 2,2) case these 
estimates are known to be sharp for Laplacian eigenf unctions as shown by Burq. 
Gerard and Tzetkov J^j. 

In proving the semiclassical version for the full manifold estimates both 
Koch, Tataru and Zworski [9] and Burq, Gerard and Tzvetkov [2] used the 
assumption (Al) along with the implicit function theorem to write as 

p(x,t) = e(x,Q(t 1 -a(x,t))- (6) 

Then by using x\ as a time variable, t, they reduced the problem to studying 
the evolution equation 

(hD t - a(t,x',hD x >))u = 0. 

An approximate propagator U(t) for ([lj can be written down as a Fourier In- 
tegral Operator. By proving a decay estimate on \\U(syU(t)\\ Ll ^ Laa they were 
able to use Strichartz estimates to determine the mixed "space-time" norm. Us- 
ing the Strichartz estimate for the pair (p,p) they obtained an estimate on the 
LP norm for p = 2 ^2^P ■ Prom the localisation assumption and Duhamels prin- 
ciple they determined the L°° estimate. All other LP estimates were obtained 
by interpolation between these points and the trivial L? bound. 

We follow a similar procedure to find estimates for |M| £P (y). As the L°° 
estimate on the submanifold must be the same as over the full manifold we only 
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need to find the L 2 norm and the LP norm given by the appropriate Strichartz 
estimates. 

We cannot however use this method immediately, as we do not know whether 
the time variable t — x\ determined by ^ remains a valid co-ordinate when 
restricted to the submanifold Y. For example, t could be constant on Y. How- 
ever the localisation property comes to our aid at this point and allows us to 
prove the required estimates (or better) when t is constant on Y . This provides 
a natural division of the problem into two cases. In case one the time variable 
is constant on Y and, given the symbol factorisation, the proof of Theorem |1.7| 
follows easily from conservation of energy and localisation. In the second case, 
where time is a coordinate when restricted to Y we need to use Strichartz esti- 
mates. Although the usual form of Strichartz estimates do not fit this problem 
we are able to modify the abstract Strichartz estimates for our use. 

The usual statement of Strichartz estimates assumes L? boundedness. In this 
case our family of operators W(t) will be determined from the the full evolution 
operator by a restriction of some spatial variables and therefore is not necessarily 
L 2 bounded. However in the Keel-Tao [7] picture of Strichartz estimates which 
we will use this unitarity does not matter. We need only to have a bound from 
which to interpolate. Obviously having a different interpolation endpoint will 
somewhat change the relationship between the Strichartz pair (r,p) and n. 

As we have shown that areas where |p(x,£)| > 1/C make negligible contri- 
butions we can study p(x,£) around the points (xo,£o) where p(xo,£o) = 0. 

In Section [2] we will factorise the symbol to create an evolution equation 
and show that if d^ z p(xo,S,o) ^ the localisation condition is enough to prove 
Theorem 1 1.7[ Section [3] gives the necessary extension of the abstract Strichartz 
estimates and governing equation for the Strichartz pairs (r,p). Section [4] uses a 
Fourier Integral Operator to represent the evolution operator U(t) and obtains 
estimates for the restriction of U(t) to the submanifold. Section [5] uses the 
estimates from Section 4 with the adjusted Strichartz estimate to prove Theorem 

o 

2 Symbol Factorisation 

By assumption (Al) we have that when p(x, £) = 0, then J^- ^ for some 
i. By the implicit function theorem we can solve the equation £j = a(x, £') on 
{£ I P( X Q>C) — 0} and, on the support of Xi we have 

P(x,0 = e(x,£)(& - a(x,C)) 

where e(xo,£o) 7^ 0. Now, as u is a quasimodc, 

e(x, hD)(D Xi — <z(x, hD x i))u — O^ih). 

As e(xo,£o) / we can, locally, approximately invert e(x,hD) so now we have 
that, 

(hD Xi - a(x, hD x ,))u = L i(h). 
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We the study the associated homogeneus evolution equation 



(hD x% + a(x,hD x ,))u = 

where the xi space variable is thought of as the "time" variable. If we can 
understand the properties of the evolution operator U(t) we will then be able 
to use Duhamel's principle to obtain estimates for u. 

As we are estimating the restriction of u to a submanifold we want to study 
a restricted for of U(t) defined by 

W(t) = R Y o U{t). 

It is now important to determine whether our time variable is a "z" variable 
(ie Y is contained in a single time slice) or a "y" variable (ie Y is transverse 
to time slices). To deal with this we will split the proof of Theorem 1.7 into 
two cases. Case 1, where d^^p(xo,£,o) (the easy case) is proved below. Case 
2, d^ r p(xo, £o) = (the harder case) requires the use of abstract Strichartz 
estimates that allow for non-unitary energy bounds. 

Proof of Theorem \1.7\ in Case 1. We will prove that if d^ z p(xo, £o) 7^ the LP 
estimates for u are at least as good (and possibly better) than those given by 
Theorem jl.7| This assumption implies ^ for some i; we assume i = 1. 
We can therefore factorise the symbol as 

where z = (z\, z'). As Pu = 1,2 (hi) and e(xo, £0) / Owe can conclude that 

(hD Zl -a(x,hD y ,hD z ,))u = hf(x) (7) 

where \\f\\ L 2 — 0(1). The associated homogeneous evolution equation is 

(hD Zl -a(x 7 hD y ,hD z ,))u = 0. (8) 

Now allowing the variable Z\ to act as a time variable we can find a propagator 
U (t) that gives a solution for d8l . The solution operator U (t) will be unitary on 
L 2 . 

Using Duhamel's principle and denoting x' — (y, z') we write 

f Zl 

u(z\, x ) = U(zi)u(0, x ) + i / U(z\ — s)f(s,x)ds. (9) 







Combining ^ with the conservation of L? mass for the homogeneous problem 
we have that if u is L 2 normalised the the L 2 mass of u on the hypersurfacc 
H = {x\zi = 0} is of order one. We now use the localisation assumption along 



with semiclassical Sobolev estimates (Lemma 1 .4 1 to obtain an estimate for the 
L? on the submanifold Y. 

HvMi* < Ki/.o.sOLoom S h- 1 ^ 1 \\u(y,0,z')\\ LlL2 
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Which (apart from the hypersurface case where it is better) is the estimate we 
are looking for. □ 

Therefore without loss of generality we will, for the rest of this paper, assume 



d^ z p(xo 1 £o) = (case 2) which by (Al) from Definition 1.6 implies d^ y p(xo, fo) 7^ 
0. To prove the estimate in this case we use the same kind of symbol factorisation 
but this time y\ will act as the time variable. 

We will use a Fourier Integral Operator representation of U(t) to obtain 
L 2 — ► L 2 and L 1 — > L°° bounds for W(t)W*(s). We can then use the Strichartz 
estimates to get an estimate on 

\W{t)u\ r p dtj' r 

where p = r. However as we will be fixing some of the spatial variable at zero 
we cannot guarantee that W(t) will still be unitary. To deal with this we need 
to make an adjustment to the abstract Strichartz estimates. 



3 Extended Strichartz estimates 

Working with the Keel-Tao [7 formalism we have a family of operators W(t) 
such that 

W(t) : H -> L 2 {X) 

for some Hilbert space H and measure space X. When we apply this we will 
have H = L 2 (R n_1 ) and X = R fe_1 . Note that M™" 1 is a time slice in M 
and X = K fc_1 is a time slice in Y. The Strichartz assumptions modified to 
include a semiclassical parameter h (see Koch-Tataru-Zworski [5] and Burq- 
Gerard-Tzvetkov 0) are that, 

\\W(t)f\\ L2 <C\\f\\ H 

and 

\\W(t)W*(s)f\\ LOO(x) < h-"(h +\t- s\)-° \\f I . 
This gives a mixed norm estimate of 

(/ \\w(t)f\\ r LP dt\ 1/r <h-&m B 

where 

2 2(7 

— I =<J 

r p 

and (r,p) ^ (2, 00). 

We adjust these estimates by allowing the L 2 norm of W(t)W*(s)f to have 
a bound of a similar form to the L°° bound. 
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Proposition 3.1. Let W(t), telie a family of operators W(t) : H — > L 2 (X), 
where H is a Hilbert space and (X, dx) is a measure space. Assume that W(t) 
satisfies the estimates 

• For all t,s eR and f G L 1 (X) 

W{t)W*{s)f\\ Loa{x) < h-»-(h +\t- s\)-°°° \\f\\ L1(x) (10) 



• For all i,s,el and f S i 2 (X) 

ir(*)w*( S )/ii i2W < + 1* - s\r° 2 

then 

\\W(t)ff Lp dt) <h W"°o — 2) + 2<«oo— 2) 

/or pairs of (r,p), 2 < r < oo, 2 < p < oo smc/i £/ia£ 

2 1, , CToo 

- + "(o-co - 0-2) = 

r p 2 



\\L 2 (X) 



II H 



(11) 



(12) 



(13) 



Proof. Following Kccl-Tao [7] we will prove the bilinear form of the estimate 



((W*(s)F(s),W*(t)G(t)) dsdt 



(14) 



Converting (10 1 and (111 into bilinear forms we have the estimates 

\{W*(s)F(s), W*{t)G{t))\ < h-»- (h+\t- \\F(s)\\ L1 \\G(t)\\ L1 

and 

'n-"(.sir ( .s..ir*"Yv ( /,) : . /, "-■(/, + \t- s\)-°* \\F( S )\\ L2 \\G(t)\\ L2 

Interpolation between these estimates yields, 

| (W*(s)F(s) > W*(t)G(t)) | < h-^^(h+\t-s\)-^^ \\F(s) llLP , ^ [mLP , 
where 

at \ 2 (°~2-°-oc) . 

V 

We now use Hardy-Littlewood-Sobolev for the t and s integrations. This 
will give us the equation governing the relationship between r and p. We have 
that 

f{x)g{y) 

_ dxdy < \\f\\ Lqi \\g\\ Lq2 

\ x y\ 



for < 7 < n, and 



1 1 _ 7 
q[ q' 2 n 
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In this case we set q\ = <?2 = r' and 



2(<T 2 - (Too) . 

7 = 1" o-c 

P 



so Hardy-Littlewood-Sobolev gives us 

2 = 2(<7 2 - goo) 
r p 

Rearranging this gives 



2 1/ N CT 00 

- + - CTOO - ^2 = — 

r p 2 
as the governing equation for these Strichartz estimates. Note that when o-i — 
= and Coo = er this is just the original abstract Strichartz estimates 
governing equation 

2 2d 

- H = (7. 

r p 

Now we need to substitute the governing equation into the h index. Doing this 
and working through the algebra we get that 



|ir(7)/|^(//) 1/r < 



_ ( Moo I er oo^2- <7 2A t oc \ 

Yl V r ( ff co -tr 2 ) 2(0-00-0-2) / # 

□ 

Note that this simplifies considerably when \x\ = <J\ and ^2 — ^2, to become 

l/r 



Remark 3.2. If is of course possible to further generalise these estimates by as- 
suming L q bounds on W(t)W*(s) for some (q , q\) rather than the usual (2, 00). 

4 Approximate Propagator 

Proposition 4.1. Suppose U(t) : L 2 (R d ) -> L 2 (M d ) satisfies 

hD t U(t) + A(t)U(t) = 0, f/(0)=/rf 

where Aft) is a pseudodifferential operator such that the symbol principal symbol 
of A(t) is real and has no dependence on h. Then there exists some t > 
independent of h such that for < t < t 

U(t)u(x) = — !— ^ / / \iM t - £ 'ri- w -rib(t,x,ri,h)u(w)dwdTi + E{t)u{x) 



(2irh) d 
where 

d t (j)(t, x, n) + a t (x, d s <j>(t, x, n)) = 0, 0(0, x, n) = x ■ n 
b(t, x, n, h) G C C °°(M x T*R d x K) £(i) = 0{h°°) : S' -» 5 
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Proof. This is in fact the normal parametrix construction yielding the eikonal 
equation for the phase function. See [5] Section 10.2. □ 

Similar to the proof in case one we will use symbol factorisation to obtain 

p(x, f) = e(x, £)(£ Wl - a(x, £')), 

where £' — (£i,£ z ) an d study the evolution equatioon 

hD t — a(x, hDyi, hD z ) = 

(see Section p| . Here we are using yi as the time variable thus our coordinate 



x is now decomposed as x — (t,y',z). In the notation of Proposition 4.1 
A{t) — a(x, hDyi, hD z ) and x — (y' , z). 

As on {£ | p(xo,£) = 0}, £ Vl = a(x,£') the second fundamental form is 
given by 

d 2 a 

The non-degeneracy condition (A2) implies is a positive definite matrix, 
therefore on a small enough patch d^a (where r\ is the dual variable to x = (j/, z)) 
is also positive definite. Recall that W(t) = Ry ° U(t) so we have (for d = n— 1) 

= T^Uf JJe^^y'^-^h^y^^h^dwdv 

In what follows we will write <f>(t, (y' ,0),r]) = 4>{t,y' \rj) and for r\ € M. d under- 
stand (y',77) = ((y', 0), 77). All dashed variable are in R k ^ 1 and all undashed 
variables are in M. d = M n_1 . 

Proposition 4.2. IfW(t) is as above then it satisfies the estimates 
\\W(t)W*(s)f\\ Loa < h-^(h+\t-s\)-^ \\f I 

\\W(t)W*(s)f\\ L2 < h-^(h+\t-s\)-^ ||/|| 2 

Proof. First we get a L°° bound on the Schwartz kernel of W(t)W*(s). This 
result can be found in [9] but for convenience we repeat it here. Using the 
integral representation for U(t) and the fact that W(t)f is the restriction of 
U(t)f to Y we write W(t)W*(s)f as, 



W(t)W*{s)f = J W(t,s,y',v')f(v')dv' 

where 

W(t,s,y',v') = , n * [ e^^-y'^-^^^-^^-^BdwdndC, 

B = B(t,s,y' 7 v',w,ri,C,h) e 5(1) n C£°(R 2+6d ). To find an estimate for 
\W(t, s, y' , v') \ we will use repeated applications of the stationary phase method. 
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First we calculate the critical points in w and £ allowing us to perform the (w, £) 
integration. The phase function (j> is stationary and non-degenerate at £ = 77, 
w = d((f)(s, v', C) and so the stationary phase method implies that 

{2nh) a J Rd 

Finally we must use stationary phase again to deal with the 77 integration. From 
the initial condition on <j> in the formulation of the parametrix we can write 

00, y', v) ~ 4>{s, v', ?y) = (y' - v', 77) + (y' - v', sF(s, y', v' , 77)) 

and so defining the phase function 4> by 

4>(t, s, 77', v', 77) = <j>(t, y', 77) - (j)(s, v', 77) 

we have that 

4>{t, s, y', v', 77) = (t - s)a(0, y' , 77) + (y' - 7/, 77 + sF{s, 77', 7/, 77)) + 0(t - s) 2 . 
So the phase is stationary when 

= dr,4>= {Id + sdr,F)(y' - v') + {t - s){d v a + 0{t - s)). (15) 
When s is small, Id+ sd. q F is invertible and this implies that at a critical point 

\y' -v'\ = 0(t-s). 

The Hessian is given by 

d 2 J = stfb/ - v', F) + (t- s)(d 2 a + 0(t - s)) 

= (t- S )(d 2 a + 0(\t\ + \ S \) 

where d 2 a — d 2 a(0, 77', 77). Here we use the non-degeneracy of d 2 a to give that 
if t and s are sufficiently small then, at a critical point, 

d 2 4>=(t-s)^(y',v',t,s, V ) (16) 

where ^f(y' ,v' ,t, s,rj) is an invertible matrix, 

det*(7/',7j',7j, s,rj) > c > 

and the elements of ^(y' , v' , t, s, 77) are smooth in all variables. So for \t — s\ > 
Mh for some suitably large M we can apply the stationary phase method to 
conclude that 

\W(t,s,y',v')\<^(l+^^j ? 
<h-i{h+\t-s\yi. 
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When \t — a\ < Mh we can use trivial estimates to show that 

\W(t,a,y',v')\ < ChT d 

<h-i{h+\t-s\)-t. 

From these estimates we can obtain the necessary bounds on the L 1 — > L°° 
norm of W(t)W(s)* . We have, 

\\W{t)W(s)*\\ L1 ^ Loo < esssup|VF(i, s, y,v)\ < Chri(h + \t- s|)"l 

For the L 2 estimate we need to use the oscillations of W(t, s, y' , v') itself. First 
note that from the critical point equation ( 15 1 we have that if \]/ — v'\ > K\t — s\ 
for some suitably large K, critical points cannot occur. In this case we can 
estimate \W(t,s,y',v')\ by nonstationary phase obtaining 

/ I ' — 'I x ~ N 
\W(t,s,y',v')\<h- d (l+ 1 -^ 1 



t - s\ 



and as \y' ~ v'\ > K\t — s 

|w(Wy)i</- "" S ' - '—— S 1 1 



In view of this we split 

W(t, s, y', v 1 ) = Wi(t, s, y', v 1 ) + W 2 (t, a, y', v 1 ) 

where 

W 1 {t, s, y', v>) = C Oj^f) W(t, a, y' , v>) 

W 2 (t, a, y', v') = (l - C (j^)) W(t, s, y', v') 
and ((r) : R — > K is a smooth cut off function 

C(r) = 

We now have 



1 \r\<lK 
Irl > 2K 



W(t)W*(a) = {W{t)W*{s)) 1 + (W(t)W*(s)) 2 

where (W(t)W*(s))i is the operator with integral kernel Wi(t, s, y' , v'). Now by 
Young's inequality 

l(W(t)W*(a)) 2 f\\ L2 <h- d (l+ ]t ^j N \\f\\ L2 J N dy' 
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it therefore only remains to deal with (W (t)W* (s)) i . 
When \t — s\ > Mh we have 

where 

ip(t, s, y',v') = 4>(t, s, y', v', T)(y', v' , t, s)) 



and r](y' \v' ,t, s) is determined by (15 1 (the implicit function theo rem guaran- 
tees that given 4> v = we can solve for r\ = r)(y',v',t, s) due to (16 1 and the 
invertibility of ty(y',v')). To exploit these oscillations we square the L 2 norm 
of (W (t)W* (s)) i / and use nonstationary phase methods. We therefore need 
derivative bounds (in y') on b(t, s,y' ,v'). Lemma 4.3 gives us that bounds on 
Dy,b(t, s,y' ,v') depend only on bounds on Dpr](y' ,v' ,t,s), < \a\ < \/3\. 

Lemma 4.3. // 



is an oscillatory integral localised around a (x dependent) nondegenerate critical 
point (p(x) = (ifi(x), . . . , <Pd(x)) an d 

\DMx)\<l H (17) 
for any multi-index a and 1 < i < n then 

e lX ^ x ^a(x,Od£ - (1 + \)ie lX * {x '^ x)) b{x) 



where b{x) obeys the bounds 

\D x b{x)\<^\ 

for any multi-index (3. 

Proof. The method of stationary phase gives us the representation of 

s^Ma^^ = (i + A) ^ e 1 ^*^)) 6(2;) 
so it remains only to check the derivative bounds of b(x). We write 
b(x) = (l + X)-i [e iX ^a(x,^ 
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where <j)(x, £) = 4>(x, £) — </>(ie, </ 3 ( a; ))- From the Morse lemma it is enough to 
prove the bounds for 

m n 

feo = E& - ^)) 2 - Efe - ^)) 2 - 



In this case we have 
dxi 

so 

where 



Now let 



-2Efe-^)%^ + 2E(,,-, i )^ ) 



dxi 



dXi 



[d Xi ip(x)] = 



{d x Mx)] T v^(x,o 



(18) 



/(A, a) = / e lA, ^a(2;,£K 



da;,- ./ \ cte da;,; 



By ( 18 ) we have 



jW**)iX*l& a (x,Qdt = -[d xt <p(x)] T I < iX -«'^hW.;>u.o<«.r.O 
and integrating by parts 



so 



d.c 



■/(A, a) ^(AJd^fVja + d^a). 



The bounds on the derivatives of 6(a;) therefore follow from (17) and stationary 
phase estimates. 

□ 



To obtain derivative bounds on r?(y', v' , t, s) we differentiate (15) in y' to 
obtain 

= (J4-i,d + 0(|t| + M)) + [d y >r){y', v\ t, s)] T (t - s)V(y', v 1 ) 
where Id]~_i ^ is the (k — l,d) dimensional matrix [7dfc_ 1 |0]. Therefore 

(t - s)[d y , V (y', v', t, s)f = -(Td k - X , d + 0(\t\ + \s\))*(y', w')" 1 . (19) 
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This gives us 



\D y ,r)(y',v',t,s)\ < 



\t-a\ 



to obtain the multi-index bound we simply differentiate (19 1 leading to 

M 



\D^{y',v',t,s)\< 



1 



\t-8\ 



for any multi-index a. So by Lemma 4.3 with A = |t — s\/h and 7 = \t — s\ 1 

/ 1 \\f>\ 
\D^b(t,s,y',v')\< i^^j 

for any multi-index (3. 
Now we have 



where 



W(t, s, u', w')f{v')f{w')dv'dw' 
W(t,s,v',w')= f Wx(t,s,y', v')Wx(t,s,y',w')dy' 



(20) 



We will estimate (20 1 via non-stationary phase estimates. The phase function 
in question is 

i>{t,s,y',v') - ip(t,s,y',w'). 
From Taylor's theorem we have that 



fe-i 



V w j Mi, «, y'- »') " s, y', w')] = Y, Jjjfriv' - w ') + °(W - w '| 2 )> 

j=l y i V 3 



written in matrix form this is 



Vy [i>(t, s, y', v') - if>{t, s, y', w')] = g^j(v' - w') + 0(\v' - w'\ 2 ). 

So we study the matrix rP,t , . As 

d v <f>(t, s, y', v, r}(y', v', t, s)) = 



and 



4>(t, s, y', v', 7]) = (p(t, y', rj) - <f>(s, v' , 77) 



dy'idv'- 
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we get 



}2„< d 



ib 

3 , / — 1 

for i, j = 1 . . . k — 1. In matrix form this is 



We already have 

a 2 n 4> = (t - s)*(y', «') = (* - «)(a*a + o(|*| + 

and 

(t - 8)[dy> V (3/, v', t, s)f = -(Id k - hd + 0(\t\ + \s\))*(y', v')~\ 

so we only need an expression for [d v >r](y', v' , t, s)] T . Differentiating (151 in v' 
gives 

o = -{id k - lld + o(\t\ + \s\)) + [d v ,ri(y', v', t, *)f(t - * W, y') 

so 

{t-s)[d vl r ] {y',v',t,s)] T = -{ld k - x , d + 0{\t\ + t/)" 1 . 
Therefore 

[^tj(i/yi 1 «)] T ^^^y,M)] = ~ (wk-M^aj-^^+odti + i*!)) 



The leading term is the upper (h — 1, k — 1) block matrix of <9^a. As is 
positive definite the matrix d ,q v , is non-degenerate. Consequently 

c\v f — ■ w 1 

|V,,[V>(t J «,i/V)-V'(f,s,y / X)]| > | f _ s | ■ 
Therefore any integration by parts of (|20|) will gain a factor of 



However each integration by parts also gains a factor of 

1 



\t-a\ 

from differentiating the symbol. Overall each integration by parts gains 

h 

\v' — w'\ 
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So we have a bound on W(t, s, v' , w') of 



\w(t, s, w >)\ < h-*\t-s\-« ( i + ^ I / C ( ^ j C ( " ) dy> 



< h~ d \t 



W -v'\\,(W -w'\ 

I 

- { d-k + i) L , K -^1 j 



and 



., \ N 
\v'— v)'\ 

h 



for all N > 0. Therefore by Holder and Young 



\\(w(t)w*(s))iff L2 < h- d \ t - sr^Wh"- 1 mi, 

2 



So for \t-s\> Mh 



\\(W(t)W*(s)) 1 f\\ L 2 < h-^\t - s\-^ \f\ L2 . 

It now remains to deal with the case \t — s\ < Mh. This can be achieved by 
scaling. In this case W%(t, s, y' , v') is only supported on the region \y' — v'\ < h. 
We have that 

\W{t lSl y',v')\<h- d . 
Using Young's inequality we obtain 

\{W{t)W*{s))if\ L2 < ||/|| L2 . 

Hence 

\{W{t)W*{s)) l f\ L2 < Ch-^ih + \t - s\)- d -=^ \\f\\ L2 . 

Putting this together with the estimates we already had for (W(t)W*(s))2 we 
obtain 

\\W(t)W*(s)f\\ L2 < Ch-^ih + \t - s\)-^ \\f\\ L2 . 

As we have used one of our original spatial variables as time we have d = n — 1. 
This completes the proof. 

□ 

Remark 4.4. In these submanifold cases it is not enough to assume, as Koch- 
Tataru-Zworski [9] did in the full manifold case, that the second fundamental 
form on {£ | p(xo,^) = 0} is merely non-degenerate. This would imply that 
d 2 a is non- degenerate, however that is not enough to guarantee that the upper 
(k — l,k — 1) block matrix of d 2 a is also non-degenerate. Therefore we cannot 
prove the L 2 — > I? estimates on W(t)W*(s) if we assume only non- degeneracy. 
Note that the L 1 — ► L°° estimate does however still hold under the weaker 
assumption of non- degeneracy. 
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We can now use Strichartz estimates (Proposition 4.2 1 on W(t). We are in 
the case that \i\ = a\ and \i2 = 02, so we have 



ww(t)f\\ r LP dty r <h-^ 



when 



1 k- 1 _ n- 1 
r 2p 4 
So this gives us that when r = p 



P 



In particular for k = d — 1, 



2(fc + l) 
n- 1 ' 



2n 
n- 1' 



Hi 2 



When k = n — 2, p = 2 so this is an endpoint. When k < d — 3, 
the Strichartz estimates give us no point (p,p). 



< 2 so 



5 Completion of Proof in Case 2 



Recall that Case 2 was d^ z p(xo,£,o) = an d so by (Al) in Definition 1.6 



d£ y p(xo,(,o) ^ 0. Without loss of generality we assume d^ yi p(xo, £o) 7^ 0. 
Around the point (xq,^o) where p(xq,£,o) = we use d^ y p(xo,^o) ^ and 
the implicit function theorem to factorise p(x,£) as 

p(x, £) = e(x, £)(f yi - a(x, £')). 

So Pit = 0^2 (ft,) implies 

e(ir, hD){hD yi — a(x, /i_Dj,/ , hD z )u = 0{h). 

As e(x, /i-D) is elliptic this implies 

(hD yi - a(x, hD y >,hD z )u = hf(y 1 ,x') 

where \\f\\ L 2 (M) = 0^(1). 

Using Duhamel's principle we write 

u(l/i, x') = C/(2/!)u(0, a/) + * [ 1 Ufa - s)f(s, x')ds. 



When we restrict to the submanifold Y by setting z = 0we get 

fVl 

u{ yi ,y', 0) = W( yi )u(0, x') + i W( yi - s)f(s, x')ds. 
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As we already have the L°° estimates we are looking for a bound for the L 2 
norm and the the bound given by the Strichartz estimates where appropriate. 
Using Minkowski's inequality we have for any q 



\\»\\l^.y) ~ ( J Wiyi^oftf, d vi 



1/3 



\ 1/9 

\W(y 1 -s)f(s,x')l q Llr dy 1 j ds (21) 
where u = u{Q,x ! ). Therefore to obtain a L q bound we need to estimate 

1/9 

\W{t)u \ q Lq dt 

y' 

In the case where k = n — 1 we obtain an estimate from Strichartz, see 
Applying adjusted form of Strichartz estimates with p = -^3r 



proposition 
we have 



4.2 



II u II LP (Y) 



< hr 1 '* \\u \y +h-^ [ \\f(s,x')\\ L2 ds 

x ' JR * 



For all other k either there is no pair (p,p) given by the Strichartz estimates or 
the pair is the endpoint pair (2, 2). 

We also need to obtain the L 2 estimates. These can be obtained directly 
from the bilinear form (14 1. 

Proposition 5.1. The following submanifold estimates hold 



ii ii <r J ^ 2 ' k < n — 3 
Mvor) ~ ^- 1/4; h = n -i 

For k = n — 2 

JV^-^), p>2 

Ml H Y)~ | log l/2 (1//l)ft -l/2 ) p = 2 - 

Proof. We will determine these bounds directly from the estimates on the bi- 
linear forms. We have that if 



J f \(W*(s)F(s), W*(t)G(t))\ < h- s \\F\ Ll 

then 

\\w(t)ff LHX) dtj <h- s ' 2 



?L 2 f^\ L *13 



\H ■ 
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Therefore using the estimate determined in Proposition (4.2 1 we need to get an 
estimate on 



h — »" 



\\F(s)\\ 2 \\G(t)\\ 2 



which by Holder is the same as estimating 



dsdt 



h — r 



(h+\t\)-^-*\F\ 2 JG\\ L2L2 
Li ' x 



Using Young's inequality this reduces to estimating 

(h+\t\)-^ 

As we are on a compact manifold and the "time" variable is actually one of our 
space variables this corresponds to estimating 

n-fc f C n-k 

hT^r I (/i + -r) _2 ^d-r. 
Jo 

Pulling the h out of the denominator and making a change of variable gives 
means this is equivalent to estimating 

[ c ' h 

<h ~-h ~ -h (1 + a)~ — da 

Jo 

rC/h 

Jo 

When k < d — 3 the integral is O(l) therefore 



J \\W(t)u\\ 



L 2 ' 



^ 2 

dt \ <h r 



Substituting this into (21 1 we get that 



\M\ LHY) <h-^- [1+ I \\f(s,x')l L ^ds 



< h i~ 



When k — n — 1 we estimate 

r C/h 







(1 + a)- 1,2 da = (l + cr) 1 ^ 



C/h 




So 
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\\W(t)uf L2 dt 



1/2 



Again substituting this estimate into (21 1 gives 



HI 



For k = n — 2 



so 



C/h 



For p > 2 we estimate 



(1 + 0-)— 3- <log(l//l) 

lkll L2 (r)<log 1/2 (l/^)^ 1/2 - 
\F{s)\iS 



L"' 



dsdt 



(h+ \t- s\)Pb>^>V 
by applying Holder and then Young we have 

h -p(p,^A) I h -iP(p,^,i) h r c/h {l+(j) -E0( P ^.i) d(7 



When p > 2 the integral is 0(1) therefore 

l/v 



(/ 



\\W(t)uf p dt 
= h-( 



< f l f(p. fi r. 1 )+Vp 



which implies the estimate 



□ 



We can now estimate the other LP norms by interpolation between these 
estimates thereby arriving at the full range of estimates. This completes the 
proof of Theorem |1.7| 



Remark 5.2. As noted in Remark 4-4 the L 1 — > L°° estimate on W(t)W*(s) 
holds if we weaken condition (A2) in definition 1.6 to require the second fun- 
damental form on {£ | p(xo,£) = 0} to be non-degenerate. From this estimate 
by Young and Hardy- Littlewood-Sobolev we can still obtain some estimates for 



small k and large p. If k < 



the the full ra nge of estimates hold. For 



k > 2^ we obtain the estimates given by Theorem 1.7 if 



P> 
p> 



4fc 

71-1 

4fc 



k>^ 

h — n ~ 1 
2 
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